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Atmospheric  Effects  on  the  Speed  of  Sound 


INTRODUCTION 


For  sound  ranging,  the  speed  of  sound  acts  as  a  scale  factor  which  trans¬ 
forms  the  measured  time  differences  at  the  microphones  into  distance  to  the 
source.  Thus  an  error  of  0.1%  in  the  speed  of  sound  results  in  a  distance  er¬ 
ror  of  0.1%,  or  100  m  in  10  km.  The  errors  due  to  the  vagaries  of  the  meteoro¬ 
logical  parameters  over  such  a  path  length  may  be  larger  than  this,  but,  with 
the  proper  model,  tend  to  be  random  and  cancel  in  the  long  run.  This  is  not 
so  for  the  systematic  errors  introduced  by  an  inaccurate  speed  of  sound  scale 
factor.  Because  of  the  systematic  and  sensitive  nature  of  sound  speed  errors, 
it  is  important  that  this  quanity  be  known  accurately  as  a  function  of  meteoro¬ 
logical  parameters. 


poses  is 


The  formula  used  to  determine  the  speed  of  sound  for  sound-ranging  pur- 

,1 


c  «=  20.06  (Tg)*'‘,  (1) 

where  c  is  the  sound  speed  and  T  is  the  absolute  sonic  temperature,  defined 

s 


3t  +  t 
v 


+  273,2  , 


where  tv  is  the  virtual  temperature  and  t  is  the  dry-bull  temperature,  both 
measured  in  °C.  The  goal  of  this  investigation  is  to  ascertain  the  validity 
of  Eq.  (1)  and  Eq.  (2)  over  the  ranges  of  temperature,  humidity,  and  pressures 
likely  to  be  encountered  in  sound  ranging. 

For  an  ideal  gas,  the  sound  speed  is  a  function  of  the  temperature  and 
composition  of  the  gas.  For  a  real  gas,  t:»e  equation  of  state  introduces  a 
correction  which  is  a  function  of  temperature,  pressure,  and  composition.  In 
addition,  the  effective  specific  heats  of  a  polyatomic  gas  depend  on  the  fre- 
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quency  of  the  sound  wave  because  of  vibrational  relaxation  phenomena.  After 
a  review  of  the  above  theory,  a  comparison  of  the  theory  with  the  present 
formula  and  with  laboratory  determinations  of  the  speed  of  sound  in  air  will 
be  made.  This  will  be  followed  by  considering  other  effects  peculiar  to  the 
atmosphere  and  a  summary  with  recommendations. 


IDEAL  -  GAS  SPEED  OF  SOUND 

The  linear  theory  (small-signal  theory)  of  sound  propagation  in  fluids 


yields 


(*■).  • 


where  p  is  the  pressure,  o  is  the  density,  and  the  subscript  S  denotes  that 
the  derivative  is  to  be  taken  with  the  entropy  held  constant  (adiabatie  pro¬ 
pagation).  Standard  thermodynamic  manipulation  yields 


c  vk  * 

where  is  the  molar  heat  capacity  at  constant  pressure. 


(w  )n  ’ 


Cv  is  the  molar  heat  capacity  at  constant  volume. 


C.  ■  T 
v 


(”)  v  ' 


and  k  is  the  isothermal  compressibility, 


1  (w\ 

~ v  W  T  ’ 


In  the  above  and  succeeding  discussion,  the  extensive  thermodynamic  variables 
such  as  entropy  S  and  volume  V  are  per  mole , 

Equation  (4)  is  correct  for  any  fluid.  For  an  ideal  gas, 

pV  »  RT,  (8) 

where  R  is  the  universal  gas  constant.  Use  of  Eq.  (7)  and  Eq.  (8)  and  the 
definition  of  density  yields 
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where  y  is  the  ratio  of  specific  heats,  and  M  ia  the  molar  mass  (mole¬ 

cular  weight).  Further,  since  for  an  ideal  gas 


C  -  C  -  R  , 
P  v 


(10) 


Eq.  (9)  becomes 


■(‘•ys 


(11) 


Thus,  for  the  ideal-gas  approximation,  and  M  are  the  only  two  parameters 
characteristic  of  the  gas  required  to  determine  the  speed  of  sound. 


For  a  mixture  of  ideal  gases,  the  values  of  Cy  and  M  are 

C  -  I  x.  C. 
v  i  i 


M  -  ZxiMi 


(12) 

(13) 


where  x^,  C^,  and  are  the  mole  fraction  (fractional  number  of  molecules,) 
the  molar  heat  at  constant  volume  and  the  molar  mass  of  the  ith  component  res¬ 
pectively.  Table  1  lists  these  quantities  for  the  primary  components  of  air. 

The  molar  heats  are  assumed  to  be  composed  of  translational  and  rotational 
energies  where  equipartition  holds.  Since  this  does  not  include  the  non-equipar- 
tional  vibrational  energies,  these  molar  heats  are  the  effective  ones  for  high 
frequencies  where  vibrational  energy  does  not  reach  equlibrium  during  the  pas¬ 
sage  of  the  sound  wave.  Therefore,  the  superscript  «  is  used.  The  molar  mas- 

12 

ses  are  based  on  the  unified  atomic  mass  scale  (  C  ■  12). 

Thus,  in  the  ideal-gas  approximation,  the  effect  of  water  vapor  on  the 
speed  of  sound  is  two- fold:  (1)  an  increase  due  to  decreasing  the  mean  molar 
mass,  and  (2)  a  decrease  due  to  increasing  the  mean  molar  heat  capacity.  Since 
the  virtual  temperature  is  defined  as  that  temperature  required  for  dry  air  to 
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yr~' 
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Table  1.  Mole  fractions,  non-vlbrational  molar  heat  capacities,  and  mo 
lar  masses  of  primary  components  of  air. 


Component 


Mol  fraction 


w  , 

c  /R 


Nitrogen 

0.7809 

2.5000 

28.0134 

Oxygen 

0.2095 

2.5000 

31.9988 

Argon 

0.0093 

1.5000 

39.948 

Carbon  Dioxide 

0.0003 

2.5000 

44.0108 

Dry  Air 

2.4907 

28.9641 

Water  vapor 

X 

3.0000 

18.0152 

have  the  ideal-gas  density  of  moist  air  at  a  particular  pressure,  the  absolute 
virtual  temperature  is 

TV  -  T  JT  ■  <u 

where  Mg  is  the  molar  mass  of  dry  air.  Thus  the  increase  in  sound  speed  (ef¬ 
fect  1)  is  characterized  by  using  the  virtual  temperature  rather  than  the  dry- 
bulb  temperature. 

The  net  effect,  including  effect  (2),  is  presently  estimated  at  75Z  of  the 

density  effect,  which  yields  Eq.  (2).  This  weighted  sonic  temperature  is 

2  3  4 

based  on  an  expression  due  to  Gutenberg  ’  , 

c  -  ca  (1  +  0.14x) ,  (15 

where  cfl  is  the  sound  speed  in  dry  air  and  x  is  the  mole  fraction  of  water  va- 


To  determine  the  constant  in  Eq.  (1),  one  must  use  the  gas  constant.  The 
value  given  in  the  1973  adjustment  to  the  fundamental  constants^,  which  in¬ 
cludes  changes  due  to  the  unified  atomic  mass  scale  and  the  1954  redefinition 


5 


of  the  ice  point  as  273.150  K,  is  shown  in  Table  2;  also  shown  are  more  re¬ 
cent  results.  We  shall  use  the  1973  value,  8314.41  ±  .25  JK  1  kmol  \ 

Furthermore,  we  will  assume  a  ±  .001  expected  error  for  the  molar  mass  due  to  a 

00 

±  50ppm  deviation  in  content  and  the  estimated  error  for  Cv  /R  of  ±  .0003. 
This  yields 

c°  =  (20.0577  ±  .0006) (t  +  273. 15)1/2  ,  (16) 

where  the  subscript  indicates  non-vibrational  specific  heats  and  the  super¬ 
script  indicates  ideal  gas,  and  t  is  the  temperature  in  °C. 

Table  2.  Recent  determinations  of  the  gas  constant. 


Reference 


Value  (J  K  1  kmol-1) 


Method 


8314.41  ±  .26 
8314.33  ±  .44 
8314.7  ±  .5 

8315.59  ±  .18 
8314.79  ±  .35 


Mean  of  pre-1962  data 
Density  (1964) 

Density  (1965) 

Sound  speed  (1974) 
Sound  speed  (1976) 


a.  Cohen  and  Taylor,  J.  Phys.  Chem.  Ref.  Data  2,663  (1973). 

b.  Rossini,  Pure  Appl.  Chem.  9,  453  (1964). 

c.  Din,  J.  Chem.  Soc.  London  1965,  829. 

d.  Quinn,  Chandler,  and  Colclough,  Nature  250,  218  (1974). 

e.  Gammon,  J.  Chem.  Phys.  64,  2556  (1976). 


Equation  (16)  is  for  dry  air.  The  addition  of  water  vapor  changes  C^ 
and  M  to 

c”/R  =  2.4907  (1  +  0. 2045x) ,  (17) 

M  =  28.9641  (l-0.3780x).  (18) 

or 
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T  -  t  1  +  0-1A59X _ 

s  "  (1  +  0.2045x)(l-0.3780x)  * 

where  Tg  is  the  effective  "sonic"  temperature.  The  binomial  expansion  of 


Eq.  (19)  yields 


c  -  c  (1  +  0.1597x)  , 
a 


which  is  larger  than  Gutenberg's  result,  and  represents  84%  of  the  density 
effect  given  by  the  virtual  temperature. 

The  amount  of  water  vapor  possible  in  the  atmosphere  is  normally  small 
so  that  the  linear  expansion  given  in  Eq.  (20)  is  probably  sufficient.  This 
is  shown  in  Figure  1,  where  the  square  root  of  the  coefficient  in  Eq.  (19)  is 
plotted  vs.  the  mole  fraction  of  water  vapor.  Plots  of  this  nature  are  used 
throughout  this  work.  Let  us  define 

c  -  co(l  +  6),  (21) 

such  that  6  is  the  fractional  difference  between  the  sound  speed  and  some 
reference  sound  speed,  6  *  (c  -  co)/cq.  Thus  Fig.  1  is  a  plot  of  6  vs  x. 


REAL-GAS  CORRECTIONS 


Air  at  one  atmosphere  of  pressure  deviates  slightly  from  an  ideal  gas. 
This  deviation  results  in  a  correction  to  the  compressibility  and  both  molar 
heat  capacities.  In  addition,  the  deviation  depends  on  the  amount  of  water 
vapor  present.  The  equation  of  state  may  be  expressed  as  a  virial  expansion. 


2l  . 

RT 


+ 


(22) 


where  the  first  coefficient  (unity)  is  the  ideal-gas  equation  of  state,  and 
B  is  the  second  virial  coefficient.  For  the  small  deviations  in  the  atmo¬ 
sphere,  the  expansion  may  be  terminated  with  the  second  virial  coefficient, 
which  itself  is  a  function  of  temperature. 

Using  Eq.  (21),  the  effect  on  the  speed  of  sound  may  be  determined  as 

'  *  '°  {1  +  rt[b  +  t2£i]}  ’  <23) 

where  the  superscripts  0  denote  zero-pressure  (ideal-gas)  values.  For  a  mix¬ 
ture  of  two  gases,  the  second  virial  coefficient  is  represented  by  7 

B  -  (1-x)2  Bj.  +  2x(l-x)  B12  +  x2  B2  ,  (24) 

where  x  is  the  mole  fraction  of  the  second  gas,  B^  and  B2  are  the  second  virial 
coefficients  of  the  first  and  second  gases,  and  B^  Ts  an  empirically  deter¬ 
mined  interaction  virial  coefficient. 

g 

The  International  Joint  Committee  on  Psychrometric  Data  has  recommended 
the  following  second  virial  coefficients  for  dry  air  and  water  vapor: 


Bx  -  40.70  -  13116/T  -  12  x  107/T3  ,  (25) 

B, -  29.53  -  .00669  T[  1-exp (-4416. 5/T)  J-17546/T  -  95300/T2  -  8.515  x  10?/T3  , 


B2  =  33.97-<55306/T)(1072000/T  )  ,  (27) 

where  the  units  of  B  are  cm^/gmol  and  T  =  t  +  273.16.  The  differentiation  of 

2  2  2 

Eq.  (25)- (27)  yields  the  required  values  TdB/dT  and  T  d  B/dT  .  If  the  new 


variables 


ao  -  Bi  +  <*°-1)T  dr  + 


-[ 


o  .  .2 


aL  =  2  B12  +  (y  -1)1 


n  7 

a2  -  B2  +  (V°-1)I  + 


o  ,,2 


o  ,,2 


„2  A12 


■*o]- 


-  2&1  +  aQ  , 


are  defined,  then 


c  *  c 

00  oo 


•  D + « + v + v2>] 


Table  3  lists  the  values  for  B^  B12>  and  B2  as  functions  of  temperature, 

while  Table  4  lists  the  values  of  aQ,  a^  and  a2.  Division  of  the  a’s  by 

3  -1  -1 

RT,  where  R  *  82.056  cm  atm  gmol  K  ,  yields 

c«  “  +  P*(bo  +  blX+  b2x2)]  »  (32) 

where  the  b's  may  be  expressed  by  the  following  polynomial  regressions, 

b  ,  =  0.445/T  -  76.7/T2  -  C950/T3  ,  (33) 

o  RT 

bl  “  RT  ”  “0.481/T  ,  (34) 

b2  -  -  -  -i~-~  exp  (1.91  +  960/T  +  1.77  x  105/T2),  (35) 


and  p  is  the  pressure  expressed  in  atmospheres. 

A  plot  from  -60°C  to  60°C  of  bQ  +  b^x  +  b2x  is  shown  in  Fig.  2.  The 

various  curves  at  the  higher  temperatures  are  for  0%,  25%,  50%,  75%,  and  10 0% 

9 

relative  humidity  (R.H.).  The  mole  fractions  were  calculated  from  . 


Table  3.  Second  virial  coefficients  for  moist  air. 
Taken  from  Goff'. 


T 

B1 

B. 

L2 

B 

2 

(°C) 

3 

(cm  /gmol) 

3 

(cm  /gmol) 

3 

(cm  /gmol) 

-90 

-50.4 

+ 

3.9 

-84.2 

+ 

8.6 

-80 

-43.9 

+ 

3.3 

-77.0 

+ 

8.2 

-70 

-38.2 

+ 

2.9 

-70.7 

+ 

7.8 

-60 

-33.2 

± 

2.5 

-65.1 

± 

7.5 

-50 

-28.9 

± 

2.2 

-60.2 

± 

7.2 

-40 

-25.0 

± 

1.9 

-55.8 

± 

7.0 

-30 

-21.6 

± 

1.7 

-51.8 

± 

6.7 

-20 

-18.5 

± 

1.5 

-48.2 

± 

6.5 

-10 

-15.7 

± 

1.3 

-45.0 

± 

6.3 

-2300  ± 

1900 

0 

-13.2 

± 

1.2 

-42.0 

± 

6.1 

-1830  ± 

800 

10 

-10.9 

± 

1.1 

-39.3 

± 

6.0 

-1510  ± 

400 

20 

-  8.8 

± 

1.0 

-36.8 

± 

5.8 

-1260  ± 

210 

30 

-  6.9 

± 

0.9 

-34.5 

± 

5.7 

-1074  ± 

116 

40 

-  5.1 

± 

0.8 

-32.3 

± 

5.5 

-  924  ± 

66 

50 

-  3.4 

± 

0.7 

-30.4 

± 

5.4 

-  803  ± 

40 

60 

-  1.9 

± 

0.6 

-28.5 

± 

5.3 

-  705  ± 

25 

70 

-  0.5 

± 

0.6 

-26.8 

± 

5.2 

-  625  ± 

16 

80 

+  0.8 

± 

0.5 

-25.2 

± 

5.1 

-  558  ± 

10 

90 

+  2.1 

+ 

0.5 

-23.7 

± 

5.0 

-  501  ± 

7 

11 


Table  4.  Coefficients  for  the  real-gas  correction  to  the  speed  of  sound 
in  moist  air.  c  *  c°[l  +  (p/RT)(aQ  +  a^x  +  a2X^)]. 


T 

a 

( 

0 

a 

i 

a 

2 

(°C) 

2 

(cm  /gmol) 

3 

(cm  /gmol) 

3 

(cm  /gmol) 

-90 

-19.9 

± 

11.7 

-41 

32 

-80 

-15.7 

± 

9.5 

-41 

+ 

28 

-70 

-12.1 

± 

7.9 

-42 

+ 

26 

-60 

-  9.1 

± 

5.4 

-40 

± 

22 

-50 

-  6.6 

± 

5.4 

-40 

± 

20 

-40 

-  3.9 

± 

4.8 

-41 

± 

20 

i 

U) 

o 

-  1.8 

± 

4.2 

-40 

+ 

19 

-20 

0.2 

± 

3.7 

-40 

+ 

16 

-10 

2.1 

± 

3.2 

-40 

± 

16 

-3400  ± 

2200 

0 

3.7 

± 

3.1 

-40 

15 

-2450  + 

1060 

10 

5.0 

± 

2.7 

-38 

± 

15 

-1850  ± 

530 

20 

6.4 

± 

2.5 

-40 

± 

15 

-1410  ± 

290 

30 

7.6 

± 

2.2 

-38 

± 

14 

-1106  ± 

158 

40 

9.0 

± 

2.1 

-38 

± 

14 

-  887  ± 

91 

50 

9.9 

± 

2.0 

-37 

± 

12 

-  721  ± 

55 

60 

11.0 

± 

1.7 

-39 

± 

12 

-  595  ± 

35 

70 

11.9 

± 

1.6 

-39 

± 

12 

-  500  ± 

24 

80 

12.7 

± 

1.6 

-37 

± 

12 

-  426  ± 

18 

90 

13.9 

± 

1.5 

-39 

± 

12 

-  365  ± 

15 

12 
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lo®10  XS 


(36) 


-  20.5318  -  2939/T  -  A. 922  log10  T, 

where  x  is  the  saturated  mole  fraction.  The  error  bars  at  -50°C  and  50°C 
s 

indicate  the  uncertainty  in  the  equation  of  state  correction.  Even  though 
the  equation  of  state  correction  has  a  negative  term  that  depends  on  the 
square  of  the  humidity,  the  correction  is  small  compared  to  the  density-spe¬ 
cific  heat  correction  (Eq.  (19))  and  the  effect  on  sound  speed  is  still  al¬ 
most  linear  for  a  real  gas.  This  is  shown  in  Fig.  3,  where  6  is  plotted  vs 
relative  humidity  for  various  temperatures.  The  dotted  line  represents  the 
ideal-gas,  dry-air  speed. 

Since  the  effect  of  humidity  is  so  nearly  linear,  one  may  use 

c  =  ca  (1  +  6x)  .  (37) 

A  plot  of  6  vs  T  is  shown  in  Fig.  A.  Here  6  is  determined  by  (c(xs)-c(o))/ 

c  x  .  It  is  noted  that  €  for  the  real  gas  is  different  from  that  of  the  ideal 
os 

gas,  represented  by  the  dotted  line.  At  low  temperatures,  the  equation  of 
state  correction  reduces  6  due  to  the  large  value  of  a2>  whereas  for  high 
temperatures  a2  is  smaller  and  the  non-linear  portion  of  Eq.  (19)  increases 
6.  It  may  be  noted  that  Gutenberg's  6  ( . 1A)  is  off  scale. 


1A 


DISPERSION  CORRECTION 

There  are  aeveral  Mechanises  which  cause  the  propagation  of  sound  in  gas- 
es  to  deviate  from  purely  adiabatic  propagation.  Since  these  mechanisms  cause 
an  increase  in  entropy,  they  lead  to  absorption  of  the  wave.  Also,  since  the 
non-adiabatic  behavior  effects  the  compressibility,  there  is  dispersion,  or  a 
frequency-dependent  speed  of  sound.  Several  of  the  known  mechanisms  are  un¬ 
important  for  the  frequencies  used  in  sound  ranging,  however,  relaxation  of 
the  vibrational  specific  heats  Is  important. 

The  so-called  classical  dispersion,  due  to  viscosity  and  heat  conduction, 
includes,  for  air,  rotational  relaxation  also.  This  dispersion  has  been  mea¬ 
sured  by  Greenspan^.  For  frequencies  below  several  hundred  MHz,  the  disper¬ 
sion  depends  on  the  square  of  the  frequency  and  introduces  a  positive  correc¬ 
tion  term  which  is  smaller  than  the  uncertainty  in  the  gas  constant  for  fre¬ 
quencies  less  than  5  MHz.  Therefore,  this  effect  is  negligible  for  sound  rang 
ing. 

The  effect  of  heat  radiation  is  to  tend  to  make  the  propagation  closer 
to  Isothermal  (a  negative  correction).  The  important  parameter  is  q/u>,  where 
q  is  the  constant  used  for  Newton's  law  of  cooling  to  describe  the  radiation 
transfer  and  u>  is  the  angular  frequency  of  the  wave.  When  q/u>  is  small  the 
propagation  is  close  to  adiabatic  and  when  it  is  large,  the  propagation  is 
close  to  isothermal.  For  constant  q,  it  is  clear  that  heat  radiation  is  im¬ 
portant  at  low  frequencies.  However,  P.  W.  Smith,  Jr.^  has  shown  that  q  is 
not  independent  of  u  and,  further,  that  the  maximum  possible  value  of  q/w  is 
2.2  X  10  ^  (q  has  dimensions  of  sec  ^).  Since  this  leads  to  a  change  in  speed 
of  about  1  part  in  10^  at  the  most,  heat  radiation  is  negligible  for  all 
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frequencies. 

Vibrational  relaxation  is  a  different  matter.  Since  a  large  number  of 
collisions  are  normally  required  to  establish  vibrational  energy  equlibrium, 
high-frequency  sound  will  not  allow  the  vibration*!  energy  to  follow  the  tem¬ 
perature  variation  of  the  sound  wave,  whereas  low-frequency  sound  will  allow 
the  establishment  of  equlibrium.  Therefore  the  effective  specific  heats  are 
functions  of  frequency  and  thus  so  is  the  sound  speed.  At  very  high  frequen¬ 
cies  the  vibrational  specific  heat  is  "frozen"  out,  leading  to  the  sound  speed 
previously  calculated  c^.  In  the  other  extreme,  for  very  low  frequencies  the 
sound  speed  is  cq  where  the  vibrational  specific  heats  participate.  The  se¬ 
paration  of  frequency  regimes  is  characterized  by  the  relaxation  frequency 
fQ  of  the  process. 

In  air  this  phenomenon  is  complicated  by  the  presence  of  several  modes 

of  vibration  which  are  coupled  by  vibrational  quantum  exchange.  The  reaction 

rates  for  the  various  collision  partners  and  a  normal-mode  solution  are  re- 

12 

quired  to  determine  the  relaxation  frequencies  .  Although  this  has  been  done 

j  3  14 

for  moist  air  at  room  temperature  ’  (five  modes  of  vibration  and  24  reaction 
rates),  there  are  uncertainties  in  some  of  the  reaction  rates  which  effect  the 
dispersion  for  low  humidities  (and  thus  dispersion  in  the  sound  ranging  fre¬ 
quency  range).  In  addition,  the  temperature  dependence  of  some  of  the  more 

o  1 4 

important  reaction  rates  is  not  well  known.  Nevertheless,  20  C  calculations 

show  that  there  are  three  important  normal  modes  whose  strengths  correspond 

approximately  to  the  Nj,  O21  and  CO2  vibrational  specific  heats.  Furthermore, 

9 

it  is  the  opinion  of  Sutherland  that,  for  high  humidities,  the  1^0  vibrational 
specific  heat  relaxes  with  the  O2  mode,  rather  than  separately.  (  They  are 


1 


1 
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connected  through  the  near-resonance  transfer  of  vibrational  energy.)  Thus, 
to  a  first  approximation,  one  may  treat  the  case  of  moist  air  as  composed  of 
three  vibrational  specific  heats  with  separate  relaxation  frequencies  more- 
or-less  determined  by  relaxation,  Oj-I^O  relaxation,  and  CC^-i^O  re¬ 

laxation. 

14 

Under  these  assumptions  , 


where 


l+(f/ft)2 


Ai  "  ki/k® 


(38) 


(39) 


and  *  adiabatic  compressibility  of  the  ith  mode,  *  high-frequency  adia- 
batio  compressibility,  f^  *  relaxation  frequency  of  the  ith  mode,  and  the  sum 
is  over  1  ■  N2,  2  -  °2+H20’  and  3  *  COj  modes.  In  addition,  Ek^  «*  k^k^,  the 
difference  between  the  low-frequency  and  high-frequency  adiabatic  compressibi¬ 
lities.  Thus 


EAi“  IT'1 


C'R 


C  (C  +  C') 
v  p 


(40) 


where  C*  is  the  total  vibrational  molar  heat  capacity.  Under  the  assumption 

that  C'  <  <  C°  ,  this  becomes 
P 


EA.  -  -  -  - 

i  00  ^00 

c  c 

v  p 


E  C 


i  » 


(41) 


or 


RC, 


Ai  * 


(42) 


C  C 
v  p 


where  is  the  net  contribution  to  the  vibrational  molar  heat  capacity  due 
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to  the  ith  mode 


The  first  approximation  to  is  given  by  the  Einstein  specific  heat 


equation 


r  -  E  *1  ej  “3 


2  exp  (-Uj ) 

^  [l-exp(-Uj)2  ] 


where  x^  is  the  mole  fraction,  is  the  degeneracy,  =  0^/T,  and  0^  is  the 
vibrational  characteristic  temperature  (the  Plank-energy  of  the  oscillator 
expressed  in  Kelvin) .  The  sum  is  over  all  vibrational  modes  that  participate 
in  one  of  the  three  vibrational  specific  heat  modes.  The  values  for  air  are: 

1)  N£  mode 

xx  =  .7809  (1-x) 

8^  m  1 

0X  =  335  ;  K 

2)  O2  -  H?0  mode 

=  .2095  (1-x),  X2  =  x 

8X  =  g2  *  1 

ex  »  2239.1  K,  e2  «  2294.7  K 

3)  C02  mode 

r'l  =  x2  =  x3  =  0.0003  (1-x) 

8l  =  g3  “  i,  g2  =  2 

01  =  1997.4  K,  02  =  960.3  K,  ©3  =  3380.1  K. 

In  addition,  there  are  correction  terms (1)  interaction  between  ro¬ 
tation  and  vibration  (centripetal  stretch),  important  at  low  temperatures, 
and  (2)  non-harmonic  oscillator  terms,  important  at  high  temperatures.  Ac¬ 
cording  to  tables  of  values  presented  in  Hilsenrath,  et  al.^  these  correction 


terms  may  be  expressed  by  the  following  regressions: 


Nj  contribution  to  air: 


~  -  (l-xH-0.8  X  10~4  +  3.5  X  10-6T), 


(44) 


O2  contribution  to  air: 


C' 

R 


(l-x)(-1.9  X  10~4  +  2.4  X  10_6T), 


(45) 


H2O  contribution  to  air: 


|  -  x  (-27.4  X  10~4  +  78.9  X  10'6T  -  1.04  X  lO'V),  (46) 


where  C"  is  the  correction  to  the  molar  heat  capacity. 

The  various  A^'s  are  plotted  vs.  temperature  in  Fig.  5,  for  both  dry  air 
and  air  at  one  atmosphere  saturated  with  water  vapor.  The  plot  is  actually 


of  -  A^/2,  since 


Caj  (1  -  fLA±)t 


(47) 


where  cQ  is  the  law-frequency  sound  speed.  Since  the  scale  of  Fig.  5  is  the 
same  as  that  of  Fig.  2,  it  is  clear  that  the  dispersion  corrections  are  of 
the  same  order  as  the  real-gas  corrections. 


The  actual  velocity  of  sound  is  a  function  of  frequency, 
c  ■  c  /  1  -  (1/2  )£  A 


1  +  (f/f ±) 


(48) 


where  f^  are  the  relaxation  frequencies.  The  frequencies  have  been  approxi- 

9 

mated  by  Sutherland  as 
* 


fn  " 


^-2" [9  +  3.5  X  104x  exp(3.8  -  25/T  ,3)]. 


(49) 


and 


‘o-  P 


(-f-)  1/2  r2 


24  +  4.41  X  10  x 


f  5  *  ^*.1 - jl  .  (50) 

V  3.91  X  10  +  x  /J 


21 


■  -J..'  - - - - —  "fl-- 


Figure  5.  Contributions  of  Vibrational  Specific  Heats  to  the  sound 
speed  vs.  temperature. 


where  n  is  for  N2  and  o  is  for  02~H20.  The  p  (293/T)  term  is  the  density  cor- 

1/3 

rection.  Additional  temperature  corrections  are  the  exp(-25/T  )  term  for  the 

1/2 

N2“H20  relaxation  and  the  (293/T)  term  for  the  °2_H2°  relaxation.  The  terms 
independent  of  x  in  Eq.  (49)  and  (50)  are  due  to  the  C02  de-excitation  of  N2  and 

°2* 

Since  low  frequen  ies  are  of  interest  in  sound  ranging,  some  idea  of  the 
temperature  dependence  of  C02  terms  would  be  helpful.  Table  5  lists  some  re¬ 
cent  determinations  of  relaxation  times  for  various  C02  collisions.  These  data 

—1/3 

are  also  plotted  in  Fig.  6,  where  log  (1/t)  is  plotted  vs.  T  (a  Landau  - 
Teller  plot),  the  predicted  major  temperature  dependence  for  exchange  of  vibra¬ 
tional  energy**.  The  de-excitation  of  N2  by  CC>2, 


n2  +  co2  +  n2  +  co2 


(51) 


h 

l 


has  been  measured  from  368  to  476  K.  Fig.  6  indicates  that  the  temperature 

1/3 

dependence  of  this  reaction  is  exp(-87/T  ).  The  9  Hz  term  in  Sutherland’s 

report  apparently  comes  from  0.0003  C02  using  the  Henderson,  et  al.  data  point. 
If  this  be  true,  consideration  of  the  temperature  dependence  indicated  by  Fig. 

A 

6  requires  this  value  to  be  lowered. 


The  near-resonance  between  N9  and  Co9(v_)  suggests  that  the  relaxation  route 

^  ^  “  *i  y  ^  ^ 

for  N2~C02  collisions  is  +  C02  -*■  N2  +  C02  (v^)  followed  by  C02  (v^)  +  M  •*  C02 

(v2)  +  M  and  C02  (v2)  +  M  -*■  C02  +  M.  Fig.  6  indicates  that  the  over-all  rate  is 

somewhat  slower  than  C02  (v2)  +  M  C02  +  M,  which  would  be  the  controlling  rate 

for  a  fast  quantum  exchange.  Bass  and  Hottman  (Ref.  33)  have  recently  measured 

the  quantum  exchange  at  200K,  finding  it  about  two  orders  of  magnitude  faster 

*  A 

than  the  relaxation  of  C02  (v3)  to  C02  (v2)  due  to  CC>2  or  N2  collisions. 


ii  < 


*  3 
1  3 
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Table  5.  Recent  determinations  of  relaxation  times  due  to  C0„  collisions. 


Ref. 

T(K) 

1/t 

t-1/3 

ln(T/293t) 

N2*  +  C02  ^  N2  +  C02 

a 

476 

1.8  X  105 

.128 

12.6 

b 

448 

1.9  X  105 

.131 

12 

c 

368 

1.0  X  105 

.140 

11./ 

°2*  +  C02  -  °2  +  C02* 

d 

300 

1.6  X  105 

.149 

12.0 

450 

2.1  X  105 

.130 

12.7 

600 

2.5  X  105 

.119 

13.1 

co2*  +  n2  -  co2  +  n2 

e 

300 

.8  X  105 

.149 

11.4 

600 

5  X  105 

.119 

13.8 

C°2*  +  02  -»•  C02  +  02 

d 

300 

1.1  X  io5 

.149 

11.6 

450 

3.2  X  10S 

.130 

13.1 

600 

6.5  X  105 

.119 

14.1 

a. 

Henderson,  4th 

International  Congress  on 

Acoustics, 

Copenhagen  (1962) 

b.  Henderson, et  al. ,  J.  Acoust.  Soc.  Am.  45,  109  (1969) 

c.  Bauer  and  Schotter,  J.  Chem.  Phys.  51,  3261  (1969) 

d.  Bass,  J.  Chem.  Phys.  58,  4783  (1973) 

e.  Shields,  Warf,  and  Bass,  J.  Chem.  Phys.  58,  3837  (1973) 


in  (T/Z  93  r) 


100  So 


Figure  6.  Temperature  dependence  of  CC^  -  induced  relaxation. 
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(52) 


The  expression  used  here  for  f  is 

f n  *  p  *  ~  [2.3  exp(13.1-87/T1/3)+3.5  X  104x  exp(3.8-25/T1/3) ] . 

The  effect  of  C02  on  02  is  supposed  by  Bass"^  to  be  controlled  by  the 
quantum  exchange  reaction 

0*  +  C02  -v  02  +  CO*  ,  (53) 


followed  by  the  faster  (due  to  the  larger  number  of  02  and  N2  molecules)  de¬ 
excitation  of  C02  by  M  (02  or  N2) 

* 


C02  +  M  C02  +  M.  (54) 

1/3 

The  temperature-dependence  of  reaction  (53)  is  determined  to  be  exp(-25/T  ) 

from  Fig.  6.  Thus  we  will  use  (55) 


fQ  -  24  P*  ~  exp  (5.4-25/T1/3)  +  4.41  X  106 


(■W 


1/2  -4 

'  5  X  10  +  x 

x  -  - 

3.91  X  10  J  +  x 


One  is  left  with  the  C02  relaxation  frequency  to  be  determined.  Calcu- 
14  o 

lations  for  20  C  yield  the  approximately  linear  relation 


fc  -  1100  +  2.0  X  106x.  (56) 

The  first  term  is  due  to  N2  and  02  collisions,  the  measured  relaxation  times 

of  which  are  given  in  Table  5  and  displayed  in  Fig.  6,  the  temperature  depen- 

1/3 

dence  of  the  N2  relaxation  time  is  exp(-80/T  ),  while  that  of  02  is  exp 

1/3  1/3 

(-83/T  ).  The  weighted  mean  of  these  is  exp(-81/T  ).  The  H20-term  is 

due  to 

CO*  +  H20  -*■  C02  +  H20  (57) 

(or  quantum  exchange  between  C02  and  H20).  The  relaxation  time  data  for  this 

18 

reaction  have  been  surveyed  by  Taylor  and  Bitterman  .  Although  there  is  a 
large  spread  in  the  data,  the  relaxation  times  seem  to  decrease  with  increasing 
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temperature.  Taylor  and  Bitterman  suggest  as  the  best  fit  a  line  which  has 

1/3 

the  approximate  temperature  dependence  exp(27/T  ).  Thus 
fc  -  P*  [1100  exp(12.2-81/T1/3)  +  2.0  X  106X  exp(27/T1/3-4.1)l  •  (58) 

The  effect  of  dispersion  is  shown  in  Fig.  7,  where  6  -  [c(x)-c  ]/c  is  plotted 

o  o 

vs  %  R.H.  for  20°C.  For  a  frequency  of  500  KHz,  the  speed  is  the  high-fre¬ 
quency  speed  c^  (somewhat  above  c°  due  to  real-gas  effects),  whereas  for  a 
frequency  of  5  Hz,  the  speed  (except  at  very  low  humidity)  is  the  low-fre¬ 
quency  speed  of  sound  cQ,  below  cw  due  to  the  additional  specific  heat.  It 
is  seen  that  for  frequencies  in  between,  the  speed  is  intermediate.  For  in¬ 
stance,  the  500-Hz  speed  has  a  rapid  change  from  ca  to  an  intermediate  speed 
representing  O2  equlibrium  at  about  4%  R.H. ,  followed  by  a  much  less  rapid 
change  toward  cq  due  to  the  Nj  relaxation  at  about  60%  R.H. 

19 

The  T-23  microphone  has  a  frequency  response  from  about  12  Hz  to  24  Hz 
This  suggests  the  receiving  of  20  Hz  waves;  however,  higher  frequency  compo¬ 
nents  in  a  pulse  will  also  cause  a  response.  We  will  assume  that  atmospheric 
di  ;>ersion  will  separate  the  frequency  components  of  a  pulse  over  long-range 
propagation,  thus  resulting  in  20-Hz  reception. 

For  completely  dry  air,  the  relaxation  frequency  of  CO^  is  considerably 
higher  than  20  Hz  (at  least  down  to  -60°C),  the  relaxation  frequency  of  0^  is 
within  the  bounds  of  the  response  of  the  microphone,  and  the  relaxation  fre¬ 
quency  of  N2  is  well  below  the  microphone  response.  Since  completely  dry  air 
is  never  found  in  the  atmosphere,  perhaps  a  figure  of  5%  R.H.  is  reasonable 
for  "dry"  air.  In  this  case,  the  O2  relaxation  frequency  falls  above  the 
range  of  the  microphone  for  temperatures  greater  than  0°C,  and  the  ^  relax¬ 
ation  frequency  is  below  the  response  for  T  <  0°C  and  above  for  T  >  20°C.  For 
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Figure  7.  Dispersion  for  20°C,  Sound  speed  correction  vs 


higher  humidities,  the  effective  speed  is  close  to  cQ. 

The  effect  of  dispersion  expected  for  sound  ranging  at  low  humidities 
is  shown  in  Fig.  8,  where  6  is  plotted  vs.  X  R.H.  for  various  temperatures 
and  20  Hz.  It  is  clear  that  beyond  5%  R.H.  the  lines  are  straight,  with 
the  curvature  due  to  dispersion  being  largest  at  less  than  2%  R.H.  Figure  9, 
which  goe3  to  100%  R.H.,  shows  that  the  speeds  are  linear  with  humidity  beyond 


Figure  8.  Sound  speed  correction  for  20  Hz  and  various 
and  low  relative  humidities. 


temperatures 


COMPARISON  OF  MODELS 


The  theoretical  model  presented  above  may  be  summarized  as 


c  -  20.0577  (T  )J 

9 


where 


j1  ♦  p*  %  * 


X  +  b2x2)  ~Jl  di  /f1  +  (f/f 


i>2]  ]  , 


bQ  -  0.445/T  -  76.7/T2-8950/T3, 


bx  -  -0.481/T  , 


01219  ^  2 

b2  ■  -  —  j  -  exp(1.91  +  960/T  +  1.77  X  10/T  ),  (35) 

dL  -  (1-x) (-9.9  X  10~4  +  1.43  X  10"5T  -6.68  X  10~8T2  +  1.05  X  10_1°T^)  ,  (60) 

d2  -  4.9  X  10'4  -4.1  X  IQ-6!  +  1.7  X  10_V  +  3.7  X  10-11T3  +  x(1.07  X  10~3 

-3.4  X  10“6T  -2.96  X  10_8T2  +  1.65  X  10~10T3) ,  (61) 

d3  -  -1  X  10-5  +  1  X  10_7T  ,  (62) 

f.  -  p*[l. 72  2.25  X  10-2T  +  8.37  X  10_5T2  +  x(1.19  X  104  +  125T  -  .1585T2)] 

(63) 

*  —4  2  7  e  y  i v 

f,  -  P  [-5.2  +  . 133T  -  1.13  X  10  V  +  7.55  X  10  x  —  -  -  -  ■  * - 

3.91  X  10~3  +  x 

X  (T)'1/2]  ,  (64) 

f3  -  p*[460-7.12T  +  . 0318T2  +  x(1.79  X  106  -  1.29  X  109/T  +  3.86  X  10U/T2]  , 


T  -  (l+.1459x)T 

s  (l-.1546x  -  . 0773xz)  .  (66 

Equations  (60)- (65)  are  regression  fits  to  the  more  complicated  equations  in 
the  dispersion  section.  They  yield  results  well  within  the  uncertainty  for 
temperatures  -90°C  <  T  <  90°C. 


*  t  ’-l  ? 

.  f*  x 

-  t  V*.  -<44-  •  ViC  jfepfc 


_ :i.  i  -...f 


The  mole  fraction  of  water  vapor  may  be  determined  from  the  measurement 


of  the  virtual  temperature  or  relative  humidity. 

T  -T 
v 

V  »  ■ 

*  .  178T  * 


Using  Eq.  (14)  and  (18), 


(67) 


or. 


x 


ux 

a 

* 


(68) 


lOOp 

where  u  is  the  %  R.H.  and  x  is  the  saturated  mole  fraction  given  by  Eq.  (36). 

s 

The  model  is  subject  to  large  extrapolations  for  frequencies  below  about 
100  Hz  and  the  exclusion  of  viscothermal  dispersion  limits  the  high  frequency 
range  to  about  5  MHz/atm.  The  uncertainty  in  the  calculations  varies  with  tem¬ 
perature.  Table  6  lists  reasonable  uncertainties  due  to  the  various  effects. 
The  real-gas  errors  are  based  on  the  uncertainties  in  Table  4.  The  magnitude 
of  the  specific  heat  was  assumed  to  be  uncertain  by  the  amount  of  the  ^  vi¬ 
brational-rotational  correction,  varying  between  0.001  R  at  -90°C  to  0.002  R 
at  90°C.  The  uncertainties  due  to  lack  of  knowledge  of  the  relaxation  fre¬ 
quencies  were  determined  by  trial  calculations  at  5%  R.H.  The  relaxation  fre¬ 
quencies  were  considered  to  vary  by  ±  50%  at  0°C  and  30°C,  ±  100%  at  -30°C 
and  60°C,  and  ±  200%  at  -90°C,  -60°C,  and  90°C.  Zero  R.H.  uncertainties  would 
be  larger. 

A  comparison  of  this  model  with  others  is  displayed  in  Fig.  10.  Here  the 
display  is  Ac  in  m/sec  vs.  temperature.  The  two  values  of  relative  humidity, 

5%  and  95%  are  dispayed.  The  highest  curves  are  the  present  model,  Eq.  (1) 
and  (2).  It  is  seen  that  the  low-humidity  curve  is  consistently  high  and  that 
the  high-humidity  curve  is  only  within  ±  .05  m/sec  (the  dotted  lines)  between 
15°C  and  30°C,  being  about  0.2  m/sec  low  at  40°C.  The  more-or-less  parallel 


(°c) 


Figure  10-  Differences  in  sound  speed  for  various  models.  All  models 
compared  to  theory  developed  in  this  work,  with  5%  and  95%  R.H.  calcu¬ 
lations  displayed.  A:  Present  sound-ranging  method.  B:  Gutenberg’s 
formula.  C:  Regression  fit  good  below  50°C. 
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Table  6.  Uncertainties  In  the  theoretical  calculation  of  sound  speed. 
16  represents  that  due  to  the  Ideal-gas  speed,  RG  represents  that  due 
to  the  equation  of  state  correction  at  one  atmosphere,  CV  represents 
that  due  to  calculations  of  the  vibrational  specific  heat  terms,  and 
RF  is  due  to  supposed  errors  in  the  relaxation  frequencies. 


T 

ERROR 

(m/sec) 

<Uc) 

IG 

RG 

CV 

RF 

TOTAL 

-90 

0.01 

0.21 

0.01 

0.01 

0.21 

-60 

0.01 

0.11 

0.01 

0.02 

0.11 

-30 

0.01 

0.07 

0.02 

0.02 

0.07 

0 

0.01 

0.05 

0.02 

0.01 

0.05 

30 

0.01 

0.03 

0.02 

0.01 

0.04 

60 

0.01 

0.02 

0.02 

0.01 

0.04 

90 

0.01 

0.02 

0.03 

0.01 

0.04 

curves 

are  due  to  Gutenberg 

(see  Eq. 

(15)). 

The  humidity 

correction 

of  the 

present 

model  is  somewhat  better  than  Gutenberg's,  while 

the  use  of  273.2 

rather 

than  273.15  for  the 

ice  point 

causes 

the  parallel  displacement 

.  Also 

shewn  in  Fig.  10  is  a  regression  fit  to  Eq.  (59)  which  is  accurate  below 
50°C  and  is  represented  by 

c  -  20.0577  (1  +  A  +  xB)  ST  ,  (69) 

where 


A 

■»  —1.43 

X  10"4-1.34  X  10‘ 

'6T  -  1.118  X  10_7T2  +  3.03  X 

10~10T3  , 

(70) 

-4 

-5  2  -7 

3 

(71) 

B 

-  .1516 

+  5.86  X  10  T  - 

1.793  X  10  T4  +  2.00  X  10 

X  T  . 

Although  this  takes  less  than  1/10  calculator  ti.e  than  Eq.  (59),  it  can  be 
seen  to  agree  within  less  than  .01  m/sec  below  50°C.  This  regression,  however, 
does  not  agree  with  Eq.  (55)  for  very  low  humidities  (<3%  R.H.).  The  slope  of 
c  vs  x  (Eq.  71)  was  determined  by  least-squares  fitting  the  humidity  points 
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every  10%  from  5%  to  95%,  and  does  not  include  the  dispersion  "loops"  shown 
in  Fig.  8.  The  intercepts  of  these  linear  fits  (Eq.  70)  and  the  slopes  were 
fitted  with  a  cubic  in  T  between  -50°C  and  50°C.  The  deviation  beyond  50°C 
is  clear  in  Fig.  10. 

The  shift  between  the  present  model  and  Gutenberg's  suggests  that  the 
present  model  can  be  improved  by  lowering  the  effective  temperature  and  rais¬ 
ing  the  %  of  virtual  temperature  used.  The  results  of  such  a  simple  change 
are  shown  in  Fig.  11.  A  constant  0.2°  was  subtracted  from  the  "sonic"  tem¬ 
perature  and  80%  rather  than  75%  of  the  virtual  temperature  was  used  to  ob¬ 
tain  this  plot.  It  is  seen  to  be  within  ±  .05  ra/sec  of  theory  between  -40°  C 
and  30° C  (for  5%  R.H.),  with  50%  R.R.  values  being  within  this  range  beyond 
60°C.  The  80%  figure  was  chosen  because  it  allowed  the  50%  R.H.  curve  to  be 
flat.  Thus  the  value 

T  -  t  +  .  8At  +  273.0  ,  (72) 

8 

where  At  is  the  difference  between  the  virtual  temperature  and  the  dry-bulb 
temperature,  is  within  i  0.05  m/sec  of  about  90%  of  expected  sound-ranging 
weather. 

One  can,  with  small  additional  temperature  corrections,  do  better.  A 
quadratic  temperature  correction  to  straighten  the  5%  curve,  a  small  tilt  of 
the  resulting  5%  curve  and  the  use  of  82  1/2%  of  the  virtual  temperature  re¬ 
sults  in  Fig.  12.  Here 

T  -  t  -  1  X  10"3t  -  6  X  10~5t2  +  . 825At  +  273.0  .  (73) 

It  is  noted  that  this  correction  is  almost  as  good  as  Eq.  (69)  but  it  only 
involves  a  quadratic  for  A,  and  B  is  constant. 

The  effect  of  pressure  on  the  speed  of  sound  is  depicted  in  Fig.  13. 
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Here  the  50%  R.H.  curve  calculated  from  Eq.  (73)  is  plotted  for  0.7  and  1 
atm.  It  is  seen  that  even  for  0.7  atm,  the  curve  is  still  within  ±  .05  m/sec. 
The  effect  of  dispersion  errors  are  shown  in  Fig.  14,  where  the  50%  R.H.  curve 
calculated  from  Eq.  (73)  is  plotted  vs  the  sound  speed  for  0.2,  20,  and  2000 
Hz.  Again,  over  this  large  range  of  frequency,  the  approximation  falls  with¬ 
in  +  0.05  m/sec  of  theory. 

The  same  insensitivity  to  frequency  does  not  hold  at  low  humidities. 

Fig.  15  shows  that  the  calculated  speed  using  Eq.  (73)  is  about  0.1  m/sec  low 
for  2000  Hz  sound  at  10°C  and  5%  R.H.  At  that  temperature  and  humidity,  the 
O2  begins  to  relax  and  the  speeds  calculated  for  higher  temperatures  are  bet¬ 
ter.  The  dashed  line  represents  the  difference  between  the  model  and  the  0% 
R.H.  2000-Hz  speed,  being  about  0.3  m/sec  below  the  model  at  60°C  (this  is 
the  total  dispersion  difference  between  cQ  and  ck). 

Except  for  these  very  low  humidities,  it  is  observed  that  simple  cor¬ 
rections  to  the  present  model  for  sound  speed  will  agree  to  within  ±  0.05 
m/sec  of  the  more  detailed  model  over  the  range  of  temperature  and  pressure 
expected  in  sound  ranging. 


Figure  13.  Difference  between  theory  and  approximation  given  by  Eq.  (1) 
and  Eq.  (73).  Relative  humidity  -  50%,  pressure  -  1  and  0.7  atm. 


AC  (m/sec) 


■» 


Figure  14.  Difference  between  theory  and  approximation  given  by  Eq.  (1) 
and  (73).  Effec,  dispersion  curves  for  frequencies  of  0.2,  20,  and 
2000  Hz  and  re)...  _  ve  humidity  of  50%. 
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Figure  15.  Difference  between  theory  and  approximation  given  by  Eq.  (1) 
and  (73).  Effect  of  dispersion  at  5Z  R.H.  Dotted  curve  is  for  OZ  R.H. 
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EXPERIMENTAL  RESULTS 


Most  speed  of  sound  measurements  in  air  were  taken  previous  to  World 
War  II.  Open-air  measurements  are  beset  with  temperature  and  wind  difficul¬ 
ties.  Laboratory  measurements  allow  the  control  of  temperature  and  humidity. 
Data  taken  at  specific  temperatures  are  reduced  to  a  near-by  standard  tempera¬ 
ture  such  as  0°C  or  20°C  by  the  use  of  ideal-gas  temperature  dependence.  After 
the  discovery  of  dispersion,  reduction  to  low-frequency  speeds  was  sometimes 

done,  since  most  measurements  were  in  the  MHz  range. 

20 

Hardy,  Telfair,  and  Pielemeier  obtained  331.44  ±  .05  m/sec  reduced  to 
cq  at  0°C.  Their  survey  of  results  up  to  1942  yielded  a  weighted  mean  of 
331.46  ±  .05  m/sec.  The  theoretical  results  detailed  above  are  in  agreement 
with  these  values,  yielding  331.44  +  .05  m/sec.  Equation  (1)  yields  331.57 
m/sec. 

Several  investigators  have  measured  both  high-frequency  (several  MHz)  and 
low-frequency  (200  Hz-2000  Hz)  for  20°C  with  varying  amounts  of  water  vapor. 
These  results,  along  with  theoretical  curves  for  c^  and  cq  are  shown  in  Fig. 

16.  The  open  data  points  represent  data  taken  at  high  frequencies,  whereas 
the  +,  x,  and  represent  low-frequency  measurements.  The  estimated  uncer¬ 
tainties  on  individual  data  points  lie  in  the  region  of  ±  .1  m/sec.  There  is 
excellent  agreement  between  theory  and  both  sets  of  data.  Harris'  data  (+) 
were  not  reported  as  absolute,  only  relative.  In  plotting  these  data,  his  0% 
humidity  value  was  assumed  to  be  equal  to  the  c^  value  for  dry  air.  According 
to  Fig.  7,  as  the  humidity  increases,  the  data  should  fall  to  the  cQ  curve. 
Although  there  is  some  overshoot,  the  data  seem  to  do  this,  a  fact  pointed  out 
by  Harris**. 
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Figure  16.  Experimental  results  for  sound  speed  in  humid  air  at  20°C. 
High  frequency  data:  D  Reid,  Phys.  Rev.  35,  814(1930);  0  Grabau,  J. 
Acoust.  Soc.  Am.  5,  1(1933);  A  Norton,  J.  Acoust.  Soc.  Am.  7,  16(1935); 

0  Pielemeier ,  J.  Acoust.  Soc.  Am.  10,  313(1939).  Low-frequency  data: 

X  Hebb,  Phys.  Rev.  20,  89(1905)  and  Phys.  Rev.  14,  74(1919);  *Partington 
and  Shilling,  Phil.  Mi.g.  6,  920(1928);  +  Harris,  J.  Acoust.  Soc.  Am.  49, 
890(1971). 
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Plelemeier  has  also  reported  data  taken  at  30  C,  where  air  will  hold  J 

m 

s 

more  water  vapor  and  the  separation  between  cq  and  cw  is  greater.  His  data  § 

are  s'..uwn  in  Fig.  17.  The  data  were  corrected  for  a  systematic  error  he  la- 
20 

ter  reported  .  This  data  shows  good  agreement  with  the  c^  theory. 

Results  for  wide  ranges  of  temperature  are  lacking  for  air.  However, 

for  various  pure  gases,  speed  of  sound  measurements  have  been  made  over  wide 

23 

ranges  of  temperature  and  pressure.  For  example.  Gammon  has  measured  the 
speed  of  sound  in  helium  (where  there  is  no  vibrational  dispersion)  from  -175°C 
to  150°C  and  from  10  atm  to  50  atm  of  pressure.  His  data  clearly  show  that 
the  real-gas  effect  is  linear  for  pressures  in  the  order  of  one  atm.  When 
the  sound  of  speed  is  extrapolated  to  zero  pressure,  his  data  also  indicate 
that  the  deviation  of  the  ideal-gas  sound  speed  from  dependence  is  within 
about  30  ppm  over  this  range  of  temperature,  with  some  of  this  deviation  pos¬ 
sibly  due  to  inaccuracy  in  measuring  the  temperature. 

It  may  be  concluded  that  experimental  results  show  no  evidence  of  inad¬ 
equacies  in  the  theory  of  the  speed  of  sound  in  gases.  In  fact,  measurements 
of  sound  speed  are  as  accurate  determinations  of  the  gas  constant,  temperature, 
and  the  ratio  of  specific  heats  as  competing  methods.  The  only  problem  for 
air  seems  to  be  the  determination  of  vibrational  relaxation  times. 
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17.  Experimental  results  for  the  speed  of  sound  in  humid  air  at 
Pielemeier,  J.  Acoust.  Soc.  Am.  10,  313(1939). 


OTHER  ATMOSPHERIC  EFFECTS 

Sound  ranging  involves  propagation  of  signals  that  over  part  of  their 
path  have  excess  pressures  that  are  not  negligible  compared  to  the  ambient 
pressure.  In  addition,  the  propagation  includes  refraction  and/or  reflection 
from  the  ground.  Furthermore,  the  atmosphere  differs  from  a  simple  mixture  of 
gases  in  a  laboratory  bottle.  Some  of  these  phenomena  affect  the  speed  of 
the  signal. 

24 

For  a  finite-amplitude  signal,  the  speed  of  propagation  is 

c  -  cQ  +  u,  (74) 

where  u  is  the  acoustic  "particle"  velocity.  Since  u  «  p/pc,  where  p  is  the 
acoustic  overpressure,  positive-pressure  parts  of  a  signal  travel  faster  than 
the  small-signal  speed  of  sound.  Unlike  water  waves,  acoustic  waves  can't 
"break",  and  the  propagation  becomes  that  of  a  shock  wave.  The  treatment  in 
this  study  is  restricted  to  small-signal  waves:  any  finite-amplitude  correc¬ 
tions  to  the  sound  speed  must  be  treated  separately. 

The  refraction  of  sound  follows  Fermat's  principle  and  takes  the  path 
25 

of  least  time  .  The  addition  of  wind  complicates  the  determination  of  this 

path,  however,  we  assume  that  an  effective  temperature  and  wind  are  available 

for  the  solution  of  the  sound  ranging  problem.  Thus  one  is  left  with  the 

problem  of  determining  the  effective  sound  speed  from  the  effective  tempera- 

* 

ture,  which  is  the  concern  of  this  work. 

* 

Finite  amplitude  and  refraction  effects  on  the  speed  of  sound  for 
sound  ranging  purposes  are  treated  in  reference  32. 
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Due  to  both  the  source  and  recei /er  being  close  to  the  ground,  sound 

ranging  involves  direct,  reflected,  ground,  and  surface  waves  close  to  the 

air-ground  boundary.  This  type  of  propagation  has  been  investigated  by 
26 

Embleton,  et  al.  ,  with  the  result  that,  although  the  amplitude  of  the  sig¬ 
nal  is  affected,  the  speed  of  sound  is  independent  of  such  complications. 

Even  an  isothermal,  non-moving  atmosphere  has  a  vertical  density  gradient 

27 

due  to  gravity.  Some  years  ago  Bergmann  showed  that  such  a  situation  leads 
to  dispersion: 

(75) 


c  =  Co  ^  ^  * 

°  2ttX 


where 


For  an  isothermal  atmosphere,  the  density  profile  is 

p(z)  =  exp(-z/H), 


(76) 


(77) 


■  v- 


where  z  is  the  altitude,  H  =  RT/gM  (the  scale  height),  and  where  g  is  the  ac- 

2 

celeration  due  to  gravity.  This  leads  to  n  =■  (5/4) /H  ,  or  using  H  =  8  km, 

c  =  cQ(l  -  2.6  X  10_5/f2)  .  (78) 

Thus  the  dispersion  is  only  important  for  frequencies  less  than  1  Hz  (wave¬ 
lengths  on  the  order  of  the  scale  height). 

28 

Turbulence  introduces  dispersion  also.  Wenzel  has  shown  that  the  speed 

2  2 

is  reduced  by  a  factor  on  the  order  of  <(Ac/c)  >,  where  <(Ac/c)  >  is  the  mean 
square  fluctuation  in  sound  speed  over  a  region  on  the  order  of  a  wavelength 
or  less.  An  estimate  of  this  effect  is  that  for  a  value  of  AT  *>  1°C,  the 
factor  is  on  the  order  of  3  ppm,  and  for  Av  *■  1  m/sec,  the  factor  is  on  the 
order  of  10  ppm.  The  value  Av  ■  1  m/sec  is  probably  a  useful  upper  bound  since 
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this  much  wind  fluctuation  would  create  quite  a  bit  of  noise.  For  instance, 

using  Bernoulli!  principle,  the  variation  in  pressure  due  to  wind  fluctuations 

2 

should  be  on  the  order  of  Ap  *  pvAv,  or  130  dyne/cm  for  a  wind  of  v  *  10 
m/sec  with  Av  ■  lm/sec.  This  would  probably  put  the  sound  ranger  out  of  busi¬ 
ness.  Thus  turbulent  dispersion  is  neglected. 

Small  aerosol  particles  may  be  assumed  to  be  in  Boltzmann  equilibrium 

with  the  air  molecules,  thus  acting  like  huge  molecules,  increasing  the  mean 
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molecular  mass.  According  to  Junges  size  distribution  of  aerosols,  the  ma¬ 
jor  effect  should  come  from  particles  in  the  0.1  to  lym  range  (about  1000  per 

3  3 

cm  ).  Assuming  a  mean  density  of  2.5  gm/cm  ,  the  effect  of  particles  in  this 

-4 

range  is  to  increase  the  mean  molar  mass  by  about  10  amu,  less  than  the  un¬ 
certainty  in  M  due  to  CO,,  variation. 

Finally,  the  effect  of  fogs  is  quite  complicated,  but  has  been  treated 
by  the  introduction  of  the  water  droplets  as  a  third  gas  along  with  water  va¬ 
por  and  dry  air.  The  result  is  a  relaxation  process  involving  viscothermal, 

mass,  and  latent  heat  transfer  between  the  liquid  and  vapor  phases  of  water. 
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Davidson  finds  that  the  dispersion  due  to  fogs  is  only  important  below  1  Hz. 
Thus  the  atmosphere  apparently  introduces  nc  speed  of  sound  effects 


(above  1  Hz)  that  can't  be  measured  in  the  laboratory.  Refraction  and  finite- 
amplitude  effects  on  sound-ranging  effective  speeds  must,  however,  be  accoun¬ 
ted  for  separately  as  effective  temperatures  and  Mach  numbers. 


SUMMARY  AND  RECOMMENDATIONS 


A  model  for  the  speed  of  sound  in  air  which  includes  the  real-gas  effects 
of  humid  air  and  vibrational  dispersion  has  been  investigated.  Other  effects 
such  as  dispersion  due  to  rotational  relaxation,  heat  radiation,  density  gra¬ 
dients,  boundary  propagation,  turbulence,  aerosols,  and  fogs  were  considered 
and  found  to  be  unimportant  for  the  frequencies  of  interest  in  sound  ranging. 
The  uncertainty  in  predicting  the  sound  speed  is  estimated  to  vary  between 
±  0.08  m/sec  at  -40°C,  ±  .05  m/sec  at  0°C,  and  ±  0.04  m/sec  at  40°C.  Experi¬ 
mental  results  in  humid  ai  at  20°C  and  30°C  are  in  excellent  agreement  with 
the  model. 

The  present  method  of  determining  the  sound  speed  for  sound  ranging  pur¬ 
poses  differs  from  this  model  by  about  +  .2  m/sec  at  -40°C,  and  by  about  +  .2 
m/sec  (5%  R.H.)  and  -.2  m/sec  (95%  R.H.)  at  40°C.  The  present  model  is  about 
.1  m/s ec  high  at  the  standard  sound-ranging  temperature  of  10°C. 


If  the  present  formula  for  "sonic"  temperature. 


3t  +  t 
v 


+  273.2 


(2) 


is  modified  to 


Tg  -  . 825ty  +  . 174t  +  273.0  -  6  X  10~5t2  ,  (73) 

then  the  predicted  sound  speed  will  fall  within  ±  0.05  m/sec  of  the  theory 
between  -60°C  and  50°C  for  relative  humidities  greater  than  about  3%.  There 
is  no  need  for  pressure  corrections  down  to  pressures  of  0.7  atm.  And,  unless 
the  relaxation  frequencies  are  off  by  factors  of  50,  dispersion  results  will 
be  within  ±  .05  m/sec  for  humidities  greater  than  5%  R.H. 

Although  the  difference  produced  by  this  change  in  evaluating  the  "sonic" 
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temperature  amounts  to  no  more  than  0.2  m/sec  within  the  temperature  range 
from  -40°C  to  40°C,  there  are  several  reasons  for  changing  the  present  pro¬ 
cedure  : 


1)  The  errors  in  determining  sound  speed  are  systematic  rather  than  ran¬ 
dom.  An  error  of  +  .2  m/sec  will  always  predict  the  source  to  be  60m  farther 
away  at  10  km,  whereas  random  errors  tend  to  cancel  out  over  the  six  micro¬ 
phones  or  from  signal  to  signal  from  the  same  source. 

2)  The  meteorological  message  "sonic"  temperature  is  rounded  to  the 

nearest  tenth  of  a  degree,  implying  ±  .05°C,  or  ±  .03  m/sec  accuracy.  Like- 

31 

wise,  tables  of  sound  speed  in  the  sound  ranging  field  manual  are  expressed 
to  the  nearest  tenth  of  a  meter  per  sec.  In  this  case,  the  implied  accuracy 
is  not  justified. 

3)  The  change  required  is  in  the  meteorological  message,  not  the  sound 
ranging  solution.  This  will  simplify  the  implementation,  since  the  difference 
in  computational  effort  between 


and 


t 

s 


. 75t  +  . 25t 
v 


t  =  . 825t  +  .174t  -  .2  -  6  X  10~5t2 

8  V 

is  not  great  for  persons  trained  to  reduce  meteorological  data. 

It  is  the  author's  recommendation  that  the  change  in  method  of  calcula¬ 
ting  the  "sonic"  temperature  be  considered  along  with  other  improvements  in  the 
sound  ranging  system. 
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